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A COMPUTER PROGRAM TO DETERMINE 
THE EFFECT OF CHARGED-PARTICLE IRRADIATION ON 
SOLAR-CELL OUTPUT POWER 


A. F. Obenschain 
Space Power Technology Branch 
Spacecraft Technology Division 


ABSTRACT 


A major consideration in designing a solar array is the charged particle 
environment in which it will operate. The computer program described here, 
developed to calculate the effects of charged-particle flux on the power output of 
a single solar cell, will greatly reduce the lengthy and laborious hand calcula- 
tions now necessary in designing an array. The program transforms the particle 
environment of a given orbit into a 1-Mev equivalent electron flux, and degrades 
an individual solar-cell I-V characteristic to account for the effect of this flux. 
The output of the program is the value of the damage-equivalent normally inci- 
dent (deni) 1-Mev electron flux and a series of current- voltage points represent- 
ing the I-V characteristic of the degraded solar cell. 
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A COMPUTER PROGRAM TO DETERMINE 
THE EFFECT OF CHARGED-PARTICLE IRRADIATION ON 
SOLAR-CELL OUTPUT POWER 


I. INTRODUCTION 

A FORTRAN computer program has been developed to account for the rffect 
of charged-particle irradiation on the power output of a single solar cell. The 
program converts the charged-particle environment of a given orbit to a damage- 
equivalent normally incident (deni) 1-Mev electron flux, then degrades the I-V 
(current/voltage) characteristic of a single solar cell to account for the effects 
of the computed flux and prints out the degraded I-V curve. Cells of either 1 or 
10 ohm-cm nominal base resistivity may be used in the program. 

The program consists of three parts: 

• MAIN-reads input data, prints output data, and calls the subroutines as 
required. 

• PHI-calculates the equivalent 1-Mev flux encountered by a solar cell 
during its orbital lifetime as a function of the coverglass and the back- 
shielding selected. 

• DGRADE-degrades a solar cell I-V characteristic to account for the re- 
ductions in power caused by the equivalent flux. 

This program is based on a paper by R. Rasmussen of RCA/AED 1 and the 
Manned Orbiting Power System (MOPS) Power Computer Program^; B. F. Poinsett 
performed most of the actual programming. The program, written in the FOR- 
TRAN n computer language is compatible with the IBM 7094 computer. 


II. CALCULATION OF 1-MEV ELECTRON FLUX (SUBROUTINE PHI) 

Because the power output of a solar cell depends on the equivalent flux the 
cell encounters during its orbit lifetime, an accurate method of determining the 
effect of this flux is essential. Subroutine PHI converts to a deni 1-Mev electron 
flux the known charged-particle spectrum , of the four types of particles that 
damage solar cells most: trapped orbital electrons, trapped orbital protrons, 
solar-flare protons, and solar-flare alpha particles. Solar-cell irradiation 
degradation is usually defined in terms of this low-energy flux because the 
damage equivalency of 1-Mev electrons for each type and energy range of space 
particle has been defined. 1 Once calculated, the 1-Mev equivalent flux for each 
type of particle is summed to determine the total 1-Mev flux the cell experiences. 
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A flux table prepared for each type of particle contains the particle's flux 
population thiough specific energy ranges on a daily basis. At GSFC, G. Stas- 
sinopoulose of the Laboratory for Theoretical Studies usually makes flux tables 
for a particular orbit. Tables 1 through 4 list the energy ranges for each type 
of particle and give the population of each for a particular orbit. The electron 
environment is broken into 11 ranges; the other particle environments are di- 
vided into 18 ranges. The program will operate properly only when the particle 
populations are divided into the energy ranges specified in Tables 1 through 4. 


Table 1 

Electron- Flux Format 


Energy- 

Range 

Number 

-- 

ae 

(Mev) 

r 

N e 

(Electrons/cm 2 ) 

1 

0-0.25 

2.89 E13 

2 

0.25-0.50 

2.35 E12 

3 

0.50-0.75 

7.78 Ell 

4 

0.75-1.0 

2.5 Ell 

5 

1.0 -2.0 

1.06 E10 

6 

2.0 -3.0 

4.09 E8 

7 

CO 

• 

0 

1 

o 

3.89 E6 

8 

o 

• 

LO 

1 

o 

* 

2.95 E4 

9 

5.0 -6.0 

3.1 E2 

10 

<p 

© 

i 

• 

o 

1.8 E2 

11 

o 

• 

00 

1 

o 

• 

c- 

1.0 El 


The program converts the daily input flux tables into total flux population 
tables by multiplying the population of each energy range by the number of days 
in orbit. (The "number of days in orbit," an input parameter specified by the 
program user, is merely the time in orbit for which the equivalent flux is cal- 
culated.) The program calculates the equivalent 1-Mev electron flux for each 
type of particle by applying a damage factor (KD) to the population of each energy 
range; damage factors depend on the effective coverglass and backshielding 
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Table 2 


Proton- Flux Format 


Energy- 

Range 

Number 

AE 

(Mev) 

N p 

(Protons/cm 2 ) 

1 

0.0- 1.0 

2.32 Ell 

2 

1. 0-2.0 

3.56 E8 

3 

2.0-3.0 

4.1 E7 

4 

3. 0-4.0 

2.19 E7 

5 

4.0-5.0 

1.23 E7 

6 

5. 0-6.0 

6.84 E6 

7 

6. 0-7.0 

5.47 E6 

8 

7. 0-8.0 

3.0 E6 

9 

8. 0-9.0 

! 

1.191 E6 

10 

9.0-10.0 

1.23 E6 

11 

10.0-11.0 

1.09 E6 

12 

11.0-12.0 

1.09 E6 

13 

12.0-13.0 

1.09 E6 

14 

13.0-14.0 

4.1 E5 

15 

14.0-15.0 

3.12 E5 

16 

15.0-30.0 

2.46 E6 

17 

30.0-100.0 

2.32 E6 

18 

>100.0 

3.83 E5 







Table 3 


Solar- Flare- Proton- Flux Format 


Energy- 

Range 

Number 

AE 

(Mev) 

N f P 

(Protons/cm 2 ) 

1 

0.0- 1.0 

0 

2 

1. 0-2.0 

1.66 x 10 8 

3 

2.0-3.0 

4.1 x 10 8 

4 

3. 0-4.0 

2.19 x 10 7 

5 

4.0-5.0 

1.23 x 10 7 

6 

5.0- 6.0 

6.8 x 10 6 

7 

6.0-7.0 

5.4 x 10 6 

8 

7. 0-8.0 

3.0 x 10 6 

9 

8. 0-9.0 

1.91 x 10 6 

10 

9.0-10.0 

1.23 x 10 6 

11 

10.0-11.0 

1.09 x 10 6 

12 

11.0-12.0 

1.09 x 10 6 

13 

12.0-13.0 

1.09 x 10 6 

14 

13.0-14.0 

8.2 x 10 s 

15 

14.0-15.0 

6.8 x 10 5 

16 

15.0-30.0 

4.9 x 10 6 

17 

30.0-100.0 

4 6 x 10 6 

18 

>100.0 

5.4 x 10 5 
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Table 4 


Solar- Flare Alpha- Particle Flux Format 


Energy- 

Range 

Number 

AE 

(Mev) 

"fa 

(Alpha 

Particles/cm 2 ) 

1 

16-18 

1.09 x 10 s 

2 

18-20 

8.2 x 10 4 

3 

?0-?2 

8.2 x 10 4 

4 

22-23 

8.2 x 10 4 

5 

25-30 

1.09 x 10 5 

6 

30-32 

5.4 x 10 4 

7 

32-35 

5.4 x 10 4 

8 

35-40 

5.4 x 10 4 

9 

40-45 

5.4 x 10 4 

10 

45-47 

2.7 x 10 4 

11 

47-52 

3.2 x 10 4 

12 

52-57 

3.5 x 10 4 

13 

57-65 

1.36 x 10 4 

14 

60-80 

8.76 x 10 s 

15 

80-100 

5.4 x 10 4 

16 


6.1 x 10 4 

17 


1.23 x 10 3 

18 

>400 

2.7 x 10 3 
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densities (gm/cm 2 ) of the solar cell. Tables 5 and 6 list damage factors for 
orbital electrons and protons, based on data reported by Brown, Gabbe, and 
Rosenzweig.** Although previously trapped orbital protons were distinguished 
from solar-flare protons because of the relative uncertainty associated with the 
solar-flare proton model, both proton classifications will use the same damage 
factor table. Damage factors for alpha particles, from the work of Smith and 
Blue, ^ are based on the relationship: 


3.8 ^ (E) = 0 A (4E) 

This relationship states that an alpha particle of four times the energy of a pro- 
ton does 3.8 times the damage of the proton. Table 7 lists damage factors for 
alpha particles. The program user must specify an equivalent thickness of 
backshielding and eoverglass for each computation; Table 8 lists available shield- 
ing densities, and equivalent thicknesses of fused silica eoverglass and aluminum 


Table 5 

Damage Factors for Electrons (KDE) 


Energy 
Range No. 

i 

Shielding Number 

1.0 

2.0 

3.0 

4.0 

5.0 

1 

0 

0 

0 

0 

0 

2 

0.03 


0 


0 

3 

0.13 



0 

0 

4 

0.30 


0.10 

0.02 

0 

5 

1.15 

1.02 

0.75 

0.47 

0.25 

6 

2.70 

2.50 

2.05 

1.55 

1.10 

7 

4.15 

3.92 

3.38 

2.85 

2.15 

8 

5.30 

5.15 

4.60 

4.10 

3.30 

9 

6.15 

6.10 

5.70 

5.30 

4.60 

10 

7.30 


6.80 

6.50 

5.85 

11 

7.80 


7.70 

7.50 

7.0 
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Table 6 


Damage Factors for Protons (KDP) 


Energy 
Range No. 

Shielding Number 

1.0 

2.0 

3.0 

4.0 

l 

5.0 ! 

1 

0 

0 

0 

0 

0 

2 

0 

0 

0 

0 

0 

3 

0 

0 

0 

0 

0 

4 

0 

0 

0 

0 

0 

5 

3000 

0 

0 

0 

0 

6 

3700 

2000 

0 

o 

0 

7 

3500 

3100 

0 

0 

0 

8 

3100 

3100 

200 

0 

0 

9 

2800 

2800 

1400 

0 

0 

10 

2700 

2700 

2000 

0 

0 

11 

2600 

2600 

2100 

0 

0 

12 

2500 

2500 

2100 

100 

0 

13 

2500 

2500 

2100 

1000 

o i 

14 

2500 

2500 

2000 

1400 

0 

15 

2500 

2500 

2000 

1500 

100 j 

16 

2500 

2500 

' 

2000 

1800 

1500 i 

17 

2300 

2300 

2000 

2000 

2000 

18 

1000 

1000 

1000 

1 

1000 

1000 
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Table 8 


Shielding Numbers 


Shielding Density 
(gm/cm 2 ) 

0.033 

0.05 

0.1 

0.2 

0.3 

Equivalent mils of 
fused silica cover- 
glass 

6.0 

9.0 

18.0 

36.0 

54.0 

Equivalent mils 
of aluminum 

5.0 

7.5 

15.0 

30.0 

45.0 

Shielding Number 
for computer 
look-up 

1.0 

2.0 

! 

3.0 

| 

4.0 

5.0 


backshielding. If coverglass and backshielding with the correct thicknesses are 
not available for use in the program, the one nearest in thickness to the correct 
one should be used. 

Once the charged-particle population of each type of particle has been entered 
into the computer, and the coverglass and backshielding thicknesses have been 
selected, the computer can calculate the total equivalent 1-Mev electron flux 
using the relationship: 


<t> T - 4> e + <t > p + 4 > f P + 4> fa (Ref. l) 


where 

<t> T = total 1-Mev equivalent flux 
= electron 1-Mev equivalent flux 
4> v = proton 1-Mev equivalent flux 
4> tP = solar-flare proton 1-Mev equivalent flux 
4> (a = solar-flare alpha particle 1-Mev equivalent flux 


The contribution of each particle to the total flux is calculated separately 
by a summation process. For example, <£ p is calculated by multiplying the 
population of each of the 18 proton energy ranges by the appropriate damage 
factors for backshielding and coverglass, then summing each of these products. 
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The following expression demonstrates this summation process: 
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o 

Ae= i 


KDPxNP| cg +KDPxNP| bs I 


where 

KDP = proton damage factor for a given energy range 
NP = number of protons in a particular energy range 

The equivalent flux for each of the other types of particles is calculated similarly. 
(The value of NP above is the total flux population of a given energy range and 
not the input per-day population.) 

In summary, the input to subroutine PHI is the charged-particle environment 
of each of the four types of charged particles, flux damage factor tables, number 
of days in orbit, and effective solar cell coverglass and backshielding thicknesses. 
The population of each energy range of a given particle is multiplied by damage 
factors for both coverglass and backshielding to determine the 1-Mev equivalent 
electron flux. The individual fluxes are summed to calculate the total deni 1-Mev 
electron flux. 


IE. CALCULATION OF SOLAR CELL I-V CURVE IRRADIATION 

DEGRADATION (SUBROUTINE DGRADE) 

Subroutine DGRADE degrades a single solar cell I-V characteristic to 
compensate for the damage caused by an equivalent 1-Mev flux. The input to 
this subroutine is either the output of PHI or a user specified 1-Mev flux; in the 
latter case subroutine PHI is not called and the computer immediately enters 
subroutine DGRADE to degrade the I-V curve for the given flux. 

The undergraded solar-cell I-V curve is entered into the computer as a 
series of voltage and current points (V; and I. ) at room temperature and nomi- 
nal solar intensity (140 mev/cm 2 ). The computer interpolates the original 
points to arrive at 100 current and voltage pairs at 10-mv increments. These 
interpolated points will then be shifted to account for the charged-particle ir- 
radiation degradation. The following describes I-V characteristic modification 
to account for the equivalent flux. 
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A curve shape correction dependent on the equivalent 1-Mev flux effectively 
sharpens the knee of the curve: the greater the flux, the sharper the knee. This 
curve shape change is obtained by adding a negative "radiation resistance" (RR) 
correction to each voltage point of the curve: 

v. = v. + I. x RR (Reference 1) 

corrected initial 


where RR, in ohms, is defined as "volts at short-circuit current (VI SC) divided 
by the initial, undegraded value of short-circuit current (ISCO)." The program 
first calculates the value for VISC, based on the equivalent flux value and the 
solar cell base resistivity, according to the following equations: 


1 ->hm- cm 


VISC 


.005 




1.85 


for 10 13 < 0 T < 3 x 10 14 (Ref. 1) 


VISC = 0.01403 log lo 0 T - 0.1929 


for <£_ > 3 x 10 14 (Ref. 1) 


10 ohm -cm VISC = .0047 


log 


10 \2.5x 10 13 / 


1.43 


for 2.5 x 10 13 < c£ t < 2 x 10 15 (Ref. 1) 


VISC = 0. 0129 log 10 <£ t - 0.1852 


for 4 > t > 2 x 10 15 (Ref. 1) 


Figure 1 shows the effect of the curve shape correction. Curve 1 is the original 
I-V curve and curve 2 shows the sharpening of the knee of the curve caused by 
changed-particle irradiation. 
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Next the program shifts the total IV curve parallel to the current axis an 
amount A I to account for the reduction in short-circuit current caused by the 
equivalent flux. The program calculates the relative short-circuit current 
degradation factor (R I) by the following equations dependent on cell base resis- 
tivity and flux level: 


1 ohm -cm RI 


1 - 7.13 x 10" 8 - 31 [0 T - 4 x 10 12 ] 0 ' 451 


for 


4 - 10 12 < <fi r < 10 15 


(Ref. 1) 


12 


RI 


= 3.3 - 0.167 log 10 <^ T 


for 4> r > 10 15 (Ref. 1) 


10 ohm- cm RI 


1 - 2.78 x 10 -734 - 5 x 10"] 


for 5 x 10 12 < < 10 15 (Ref. 1) 

RI = 2.806 - 0. 1325 log, 0 e£ T 


for 0 T > 10 15 (Ref. 1) 


The current increment A I is calculated as ISCO (1 - RI). Curve 3 of Figure 2 is 
determined by subtracting AI from each current point of curve 2: I. = I . - Al . 
The point PZCD is not the actual degraded open-circuit voltage, but merely a 
temporary value of the "point of zero current" of the cell. 

A final characteristic correction is made to account for the reduction of 
open-circuit voltage in the cell caused by the equivalent flux. The relative open- 
circuit voltage degradation factor (RV) is calculated by the following equations: 


1 and 10 ohm- cm RV - 1 - 0.022 


logl ° (3* 10 12 ) 


1.67 


for 3 x 10 12 < <£ t < 10 14 (Reference 1) 


R v = 1.779 - 0.0588 log 10 <*> T 


for 4> r > 10 14 (Reference 1) 
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VOLTAGE 


PZCD PZCU 


Figure 2. Short-circuit current correction. 


v 


(unlike the previous degradation parameters, Rv does not depend on the solar-cell 
base resistivity.) The degraded open-circuit voltage (VOCR) is calculated as 
(PZCU x Rv ). The curve is then shifted parallel to the voltage axis an amount 
Av— defined as (PZCD- VOCR)— by subtracting Av from each voltage point (V. = 
Vj-Av). Curve 4 of Figure 3 shows the fully degraded solar-cell I-V 
curve. 

In addition to degrading the solar cell I-V characteristic, subroutine 
DGRADE calculates the percentage of short-circuit current and maximum power 
degradation. If this additional calculation is desired, the user must input the 
initial I-V curve values of maximum power point current and voltage. Because 
the initial value of short-circuit current in the cell is already part of the input 
data (the first current point of the input I-V curve), this parameter need not be 
defined again. 
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IV. PROGRAM USAGE 


The following describes the mechanics of using the program in "non- 
programmer’* language. 

The program is broken into two parts: a source deck and a data deck. 
Appendix A contains a complete source deck listing, whereas Appendix B 
shows a typical data deck setup. As the source deck (A) is merely read into 
the computer each time without modification, the following paragraphs will de- 
scribe only the contents and format of the data deck. (Refer also to 
Appendix B.) 


f 

I 

f 

I 
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Asterisk Data 


The initial data deck card is the asterisk data card. This card tells the 
computer that data is to be read next, and has the following format: 


* .... Column 1 
DATA .... Columns 7-10 


The remainder of the card is blank. 


Charged Particle Environment 


This section of the data deck contains a portion of the information required 
by subroutine PHI: the population, through the energy range, of each of the four 
types of charged particles. Flux tables are entered in order for electrons, pro- 
tons, solar-flare protons, and solar-flare alpha particles; a maximum energy 
range card precedes each table. This card gives the number of energy ranges 
of any particle population, and this number must agree with the number of ranges 
for each particle as shown in Tables 1 through 4. The range-card number is an 
integer, right-justified against the third column of the card. The remainder of 
the card is blank. Following the maximum-energy-range card, the population of 
the particle is read in a continuous format, with the energy ranges as stored 
parameters. That is, the computer will automatically equate the first number 
read with the population of the first energy range of the particle, the second 
number with the population of second energy range, and so forth until the entire 
flux table for the particle has been read. Erroneous results will be obtained if 
the tables are not entered in the correct order or the particle environment is not 
broken into the specified energy ranges. 

Each particle population card contains eight 10-column fields; the population 
of each range is right- justified against the tenth column of each field in "E" 
notation. If a given particle type such as solar-flare alpha particles is not con- 
sidered to exist for a given set of orbital conditions, the maximum range card 
and corresponding population data cards with a zero in each field must be entered 
into the computer. 
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Damage Factors 


The next portion of the data deck contains the damage-factor tables. This 
section is preceded by an NMAX card which gives the maximum number of 
shielding densities— five, at present- the program can handle; NMAX is an integer 
number, right- justified against the third column of the card, with the remainder 
of the card blank. Following the NMAX card, damage factors for electrons are 
read in a continuous format, with shielding density numbers and energy ranges 
as stored factors. The first 11 numbers read give the damage factor of each of 
the 11 electron ranges, assuming a shielding density of 1.0. The next 11 num- 
bers correspond to the energy ranges with a shielding number of 2.0. Therefore, 
the first 55 numbers entered after NMAX (5 shielding densities x 11 energy 
ranges) give the entire electron damage factor table. Id a similar manner, the 
next 90 numbers (5 shielding densities x 18 energy ranges) make up the proton 
damage factor table. Trapped orbital and solar-flare protons both use the same 
damage factors. The last 90 numbers read make-up the alpha particle damage 
factor table. 

Each card in this portion of the data deck contains eight 10-column fields; 
the damage factors are real numbers and are right- justified against the tenth 
column of each field. As the program can only use the shielding densities of 
Table 8, additional damage factor tables must be generated if other shielding 
thicknesses are required. However, satisfactory results may eften be obtained 
from use of one of the available shielding thicknesses. 

Run Information 


The last section of the data deck contains the run information. Each card 
has the following format: 

1. NRUN— rim number, columns 1-3 (this number is an integer, right- 
justified against the third column, and is used to identify a particular 
run) 

2. BS— backshielding thickness in mils, columns 4-13 (this real number is 
right- justified against column 13) 

3. CG— coverglass thickness in mils, columns 14-23 (this real number is 
right- justified against column 23) 

4. RHO— solar-cell base resistivity in ohm-cm, columns 24-33 (this real 
number, either 1.0 or 10.0 is right- justified against column 33) 
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5. PHIT— total 1-Mev equivalent flux if read as data, columns 34-43 (a 
real number in "E" notation is right-justified against column 43 — these 
columns are left blank if no flux is to be read) 

6. NPHI— signals computer that PHIT is read as input data, column 46 
(the integer 1 is punched in column 46 whenever the total 1-Mev flux is 
read into the computer rather than being calculated — subroutine PHI 
will not be called) 

7. NCELL— signals computer that a new solar cell I-V curve will be read 
in, column 49 (the integer 1 is punched in column 49; if the same solar 
cell I-V curve is to be used, column 49 is left blank) 

8. DAYS— the actual number of days the solar cell is exposed to the 
charged-particle environment — columns 50-56 (this real number is left 
justified against column 50) 

9. XPV— the solar cell I-V curve initial, undegraded maximum power- 
point voltage in volts — columns 60-69 (this real number is left- justified 
against column 60) 

10. XPI— the I-V curve's initial maximum power-point current in amperes — 
columns 70-79 (a real number left- justified against column 70) 


NOTE 


If XPV and XPI are not entered as input data, the per- 
centage degradation of maximum power and short- 
circuit current is not calculated. 


When a run-information card indicates that a new solar-cell I-V curve is to 
be entered, the next card will be a solar- cell header card. This card contains 
the following information: 


1. NTBL— the number of the I-V curve — columns 1-3 (this integer number 
is right- justified against the third column of the card) 

2. NDTS— the number of current-voltage pairs — columns 4-5 (a 2-digit 
integer number) 

3. DATE— a 6-digit date number (e.g. 020569 = February 5, 1969) columns 
11-16 
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4. TITLE— descriptive title for solar cell — columns 21-80 (any alpha- 
numeric combination may be used to describe the solar-cell I-V curve) 

The cards containing the cell I-V points follow the header card. The data are 
entered alternately— current value, voltage value, etc — starting with the short- 
circuit current. One data point requires 10 columns; thus, all 80 columns of one 
card will handle four I-V pairs. Each number is right-justified against the tenth 
column of its field. 

For additional runs, merely add more run-information cards; each card will 
have the format explained in the foregoing. (The last portion of Appendix B 
shows a typical three- run data deck. Most of the possible user options have 
been demonstrated in this example.) The last card of the data deck contains the 
integer number 999 in its first three columns; the remainder of the card is blank. 
This card signifies that all data have been entered into the computer. 


V. PROGRAM OUTPUT 

The output of the computer program can De divided into three basic parts: 
input data, 1-Mev equivalent electron flux, and the irradiation-degraded I-V 
curve. The first section, containing all pertinent input data, is printed out with 
each run to give the user a quick and accurate tabulation of the data used by the 
program. 

Input data included in the printout are: 

1. Run number 

2. Number of days in space for which the calculation is being made 

3. Backs hielding thickness in mils of aluminum 

4. Coverglass thickness in mils of fused silica. 

5. The population of each energy range of the four types of particles 

6. Damage factors for each type of particle 

7. The original I-V curve, as a series of current-voltage pairs. 

The second section of the program output contains the equivalent 1-Mev flux 
for the four types of particles and the total 1-Mev equivalent electron flux ex- 
perienced by the solar cell during a given number of days in orbit. Section three 
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is a printout of th<? solar cell's degraded I-V curve as a series of current-voltage 
pairs and the cell's degraded values of open circuit voltage, short-circuit current, 
and maximum power current and voltage. If the undegraded values of solar cell 
maximum power current and voltage are given as input data, the percentage de- 
crease in short circuit current and maximum power is also calculated and 
printed out. 

The program can also plot the input and degraded I-V curves for each run. 
These plots will be made even if the same input I-V curves is used for several 
consecutive runs. Although these plots are more of a qualitative than a quanti- 
tative nature, they may prove helpful in a "quick-look” analysis to determine the 
relative decrease in power or current caused by various flux environments. 


VI. SUMMARY 

A computer program has been described which calculates a damage-equivalent, 
normally incident, 1-Mev electron flux from a given charged particle environ- 
ment. The program will degrade a single solar-cell I-V curve to account for 
the equivalent flux. In addition, the program will permit the repetitive calcula- 
tions necessary to determine the effect of different charged-particle environ- 
ments and various coverglass and backshielding densities on the equivalent 
1-Mev electron flux. This information is essential to the adequate design of a 
solar array. 
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GLOSSARY 


• AISC— interpolated current point of I-V curve 

• B— backshielding thickness in mils 

• C— coverglass thickness in mils 

• CISC— original input data current point of I-V curve 

• DA— solar flare alpha-particle damage factor; this parameter is called KDA 
in program writeup 

• DATE— arbitrary 6-number figure, usually signifying the date on which a 
solar-cell I-V curve table was made up 

• DAYS— number of days a solar cell is exposed to a given charged particle 
environment 

• DE— Electron damage factor; parameter referred to as KDE in the program 
writeup 

• DELTAI— I-V curve short-circuit current correction value, defined as 
ISCO (I-RI) 

• DP— Proton damage factor; this parameter is called KDP in program writeup 

• E— number of electrons in an energy range; referred to as N e in the program 
writeup 

• FA— number of alpha particles in an energy range; referred to as N fa in the 
program writeup 

• FP— number of flare protons in an energy range; referred to as N f in the 
program writeup 

• I . —a current point on the solar-cell I-V curve 

• ISCO— initial, undegraded value of cell's short-circuit current 

• I-V— designation signifying current/voltage 

• KDA— damage factor for alpha particles; referred to as DA in the program 


22 



• KDE— damage factor for electrons; this parameter is called DE in the 
program 

• KDP— damage factor for protons; referred to as DP in the program 

• LEVELE— number of electron energy ranges in the flux table 

• LEVELFA— number of solar-flare alpha-particle energy ranges in the flux 
table 

• LEVELFP— number of solar-flare proton energy ranges in the flux table 

• LEVELP— number of proton ranges in the flux table 

• MOPS— Manned Orbiting Power System 

• NCELL— signifies that a new solar cell I-V curve is to be entered as input 
data 

• NE— number of electrons in an energy range; referred to as E in the 
program 

• NFA— number of flare alpha particles in an energy range; referred to as 
FA in the program 

• NFP— number of flare protons in an energy range; referred to as FP in the 
program 

• NMAX— the maximum number of shielding thicknesses the program can 
handle 

• NP— number of protcns in an energy range; referred to as P in the program 

• NPHI— an integer which signifies that an equivalent 1-Mev flux will be 
entered as input data 

• NPTS— number of input I-V curve voltage-current pairs 

• NRUNS— the number of computer runs for a given data-deck setup 

• NT BL— number of the I-V curve table 

• P— number of protons in an energy range; this parameter is called NP in 
the program writeup 
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PHIT— equivalent 1-Mev electron flux entered into computer as input data 


PWMT— temporary trial value of the maximum power 

PWR— maximum power of degraded solar-cell I-V curve 

PZCD— degraded point of zero current; a temporary parameter used in 
Subroutine DGRADE 

PZ CU— point of zero current- the initial value of open- circuit voltage 

RHO— solar-ee*l base resistivity 

RI— relative short-circuit current degradation factor 

RR— radiation resistance, defined as VISC/ISCO 

Rv— relative open-circuit voltage degradation factor 

TITLE-tile of I-V curve 

v— irradiation-degraded voltage point of I-V curve 

VCC— original input-voltage point of I-V curve 

Vj —any particular voltage point on I-V curve 

VISC— volts at zero current; an equation 

VMP— maximum power point voltage of degraded I-V curve 

VOC— interpolated voltage point of I-V curve 

VOCR— degraded value of open-circuit voltage 

XI— irradiation-degraded current point of I-V curve 

XTMP— maximum power point current of degraded I-V curve 

XPI— maximum power point current of initial, undegraded solar-cell I-V 
curve 

XPv— maximum power point voltage of undegraded solar-cell I-V curve 
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• Ai— decrease in short-circuit current caused by the equivalent flux; defined 
as ISCO-ISCR 

• A v— decrease in open-circuit voltage caused by the equivalent flux; defined 
as PZCD-VOCR 

• 4> z —equivalent 1-Mev flux of electrons 

• <£f a —equivalent 1-Mev flux of solar-flare alpha particles 

• p —equivalent 1-Mev flux of solar-flare protons 

• 4> P —equivalent 1-Mev flux of protons 

• 4> t —total 1-Mev equivalent flux 
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APPENDIX A 

* 0 ‘- 1 r 0 A / i 6 / 6 i> 

* oOl AX C i_ L L i i \ i . . i s i (\ I 1 U L' G A mJm T I P R J Lj R m • i • 

* P A 0 -j L 


X-Rj 

SAIN PPuGPA*. 

OR A;'/ I Of; L ( i i ' -■ ( 1 c : i * FP ( 1 ° ) ? F / -. { 13) * : - A ( 11 *Si * }P( IS *5 ) 9 I .'A ( 1BAJ i 
IX I ( i.i ) * V ( 10 1) 


ALL 

E P T I u 

FAT T 

rl; a;p- stl 

R T u 

AT I 

. 1 1 T 1 ■■■ • 


aF A !') 

liLPUT 

TAPE 

2 9 l O i 9 L P V ' . Li 

9 ( 

L ( I } 9 

I f I 9 L C V l~ L ; ) 


i-L- AL 

I PPL i 

TAPE 

2 9 2 L j 9 L VL Li j 

9 * 

P l I) 9 

i = i 9 L 2 VL L.R ) 



i; ( :'uT 

i' A P ii 

2 9 2 0 3 9 L i. V L i P 

9 ( 

r 'Ml) 

9 I “ 1 9 L i— " J L , I ' ) 



I !’>P j T 

T A P r 

2 9 2 0 3 * L i- V L F A 

9 ( 

r A ( I ) 

* I = i 9 LEVEL A ) 


DAT' AG 

E r A C 

TOR T 

alle: are col 

PL 

ai; i 

• 


P.F A D 

rrut 

T APE 

2 9 20 s 9 f; : 'ax 





. r i\ ' 

I P ' 1 T 

TAPE 

2*205* ((LEU 

* J) 

9 I = I 9 

LK VKLr ) *0 = 19', 

.AX ) 

NEAL 

IfiPUT 

TAPE 

2*20 5 9 ( (ufM I 

9 J ) 

*1 = 1* 

LEVFLP) f J=1*P 

- AX ) 

PE AL- 

IPP'.JT 

TAPE 

2*2059 ( ( LA ( I 

* J ) 

*1 = 1. 

LEVLFA ) * J=1 ,N 

FAX ) 

II PEAL 

I A P U I 

tapl 

2*1 * OK URL » 

b 9 

C * K 

t P ri 1 T 9 . v r , i I 9 '■ n 

i L E L L 9 b A Y u * X P V * X P I 


.. IF tCR!RR-999) **2*2 
3 IF ( N P t i I ) 21 *20 t 


'2 ■ R I T r our PUT Taps 3.2^.^'JUD. DAY s iii ♦ c 

OUTPUT TAPE 3»13 
L>C IS I - 1 9 L T. V : ‘ LP 

I r ( I - 1 1) 1 5 9 1 1} 9 i 6 

„1 L _ R I T E o'u . R' I I a P t ^ 9 i i 9 I 9 E ( i ) » i 9 r ( i ) 9 1 9 r F ; T I ) 9 I * r ^ 0 I ) 
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GO TO 14 

16 . . R I T £ 0 U T P u f 1 ' PE 3 * 2 1 1 * I * k ( i ) * 1 * f 0 \ i ) * l * T ; ■. ( i ) 

14 CVOII il'Jv. 

■RITE CUT 0^ T TAPF 3*4 
t-061 = 1 * L V L L P 
i F ( I -11 ) T * S , ' • 

V. I r 10 OUTPUT 1 APF 0> * 2 - 8 * I *{[>f(I ♦ J ) » J= 1 * UMAX ) * I * { DP ( I *J)*J*1* 

1 0-i AX ) 

GO TO 6 

9 .. r! I f - uUTPOT r*PE A * 2 ;• V 1 1 * ( w P { I * J ) * J = 1 * Mrs A X ) 

6 COT;T I. Mill 

v.nm: output tape 3*12 

0 u 7 1 = i t L E V L F A 

a R I T L Jvj 1 Pu I T APF 3 * 2 1 ^ * I * ( D A ( I * J } * J = 1 * N M A X’ ) 

7 CuK T 1 N 1 1 r . 

GO TO* 2 2 

21 ..RITE OJTPuT TAPF 3 * 2 j7 0 »NRUM$ 

22 GALL PHI (0*;.PHlr * Pt ; 1 P * Pf ! I FP ♦ Pi ; I F A ♦ PH 1 T * DO t DP * OA * F * P *FA * FP * I * J* 
1 * 1 0 A X * L £ V E L E. , L E V F L P * L E V L r A , L £ V L F P * \ IP H I ♦ D A Y S J 

G h L L Lv G R A E ( P H I T * R * X I * V • ft C F L L * X P V » X P I ) 

GO TO 11 

2 ..RITE OUTPUT TAPE 3*10 
CALL EXIT 

1 r OR, t* T ( I 3 t F 1 u • t E 1 C « - ♦ 2 I 3 * F 7 • 0 * 3 X * 2 F 1 0 • 0 ) 
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10 F 0 k . A T { 2 C n 0 A L L K 0 *•! o L u ". H L L. T L J • ) 

20 3 FORMAT (13/1 Jl iG.G ) ) 

2 04 FOKMA1 ( I3-) 

2 3 6 FORMAT (rtFir.ri) 

207 f ORMAT ( 1 3HU RUN NURFFR X3*10X*i^ii . A Y IT : : P.A >r 7.?// 

12CH RAC Ko. lit Ln if 1 l Li>-r 6 • 1 ♦ SX * 2 ?J 1 CGVFRhlAOG IN r - : I L A •• r r ; • 1 ) 

2.70 FORMAT ( 131 !v RoN NO M/.f-R 13) 

4 F ORi l A T { 5 0 X t 2 1 it 0 A M •' • r A C T OR T A A L r " S / / 3 3 X * 1 3 - * • - *. IN I N T ! U 5; / / 

1 3 6X * V H P RGG K A . ■ t # / / / / 1 4 X t 2 9 ii i 1 A UA 0 i / Alu a A F u A . L r A 1 1- A ' * o> ♦ 

2o&X*27m t>/\.'./,oL rAClGRG r^R P i\ ~/T *jM.S / 2 X * fc>f L ; ! * 0 2 X * A. it i ,L / 

3 2X * 0 HP PNG L * j, 7 X * 1 o H 0 H I L l) I No NX'.-.; !\P*li-X OH. /•-» ^ vj* * 1 7X * x or; ,-.H I c. L - i No 
a. .ij; ,i: l ap/ 2X *6i ;.‘:i ■; ihLU *4x * 3 h l .o * 7 x * o r 1 1 .0 1 7x # . r * 7x * oho .0 *7x^n'i#n* 

5 3 X * 6 H i\ U . 't b L I'? * 4 X * 2 H 1 • r - '• « / X * 3 / i 1 • 3 * 7 X * j> m J • 0 * 7 X * 2 t . A € f * 7 X * * 1 1 v » 0 ) 

2 0 S FORMA T ( 3X • I 2 * 2 X * 0 F 1 r * 3 * 3 X * I 2 • 3 >• * r F 1 ^ ) 

2 G.v FORMA T { o 3 X ? X * I 2 * j X * . F i 0 • ; ) 

12 FORr.O. T ( 1 ' X * 3 or! L-AMAC:. FACTOR 5, FOR ALPHA P A F T K Lr ■' / 2 X • 6=7 l\i i * t . Y / 

12X • 6 HP . X c ■; * l / X t 1 0 : . R.UlT".' ?:,/ 

? 2X * 6 IliV H R ; R * 4 X • ? H 1 • * ?X * 3- U 1 • F * 7X * : I ; > • A * I X • i • * 7 X * 7! 14 • O ) 

2 1 ~ ; .pis : ' A T t 7 *12* 2 X * A t 1 0 • 3 ) 

13 f 0 R : 1 a T ( v X • 4 O- : P A o T I l L P 'L • T i ™ 0 F' r NO P' Ti M r PP^GP •' M • / 

1 v '3 X » 0 H NHL 7 A / 7 i X * 0 r r LAP*! 7 • - I i-1 L P H / 1 ’ X * 6 ,: '‘0 r *”Y ♦ 2 X * ' * 

2 7 X • 6 : : j !'3.r *2X.7tiP., ;T.;R r * -X 6r.~ NOP-.. V ,; : v • r; n' 5 -: T^c, ;v,r L r-o :Yl 
32X.';:iPARTULL /i . v * ' hRFNCl •, ?P P R • IPX • r H A ' ' A * 6 X ♦ 3 ! 1 P F !< ♦ 1 1 X * 0 i A ” 

4gf *gx • : • p„ n • '-'*:,xo no. • a v * : ■ n ! . r / 
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. ; r V . 1( y ,6 r' ■ 


, 6i p. 1 : y*\ . r P . 4 X « M IS 0 . G 1 . ♦ B X . 6HMU: 

‘^•EPt4X* 


f i v « 6 f < r 

. *C '■./) 


* ; ;;.T ( i7X » 12 ^v-‘ ? "!■’ 

'.2* .»X ‘ I 

2 ♦ 4X * =: 1 7 • 2 » 8 X • I 2 » A X • E 1 ri • 3 » 8 X f 

12 ,4X*E1 

12 ) 

r.:x ; 1 ( 4 lx * 12 *4x * rr 

» ; » ' a * I 

2 * 4 X ♦ ?• ] ; . 3 * BX * I 2 ♦ 4X * 1 r. ♦ 3 ) 


- \j ! ; 

f I Pi PH Ml.: * C i 

I * P M I 2 * •- 

’ i i 1 P * P i-i I F r • * P i i I T 2 * P ; i ' f * 1 ; 2 * T;P * 0 A 

9 r * P * F A 

1 f r ? * i * J » N ’ i-f ■/ i..L . 

' *Li\ 2 L 7 

’’ » L V LC A » L i V L A P * X P» i * * \:-r. Y S ) 


id 

) ,rP! * . : 

|,r.'(.)f i ;:(0)» ‘ PDM i)A(O) 




106 IF (63 — 3 -'*0)9 9*107* I jd 
1 ; ~ 7 NBL = 4 

GO TO 111 

106 ’ i F ( E F -4 Co f 9 9 ,109*110 
109 NBS=5 

GO TO 111 
11C PRINT l'C I 

2101 FORMAT ( OX *4 1H 3ACKSH IEL0 TOO THICK - '-'ILL USE 43.0 
N3S = 3 
GO TO 111 

COV-RGCASO .INTEGER NU?«bER FOR TABLE LuOK-UP. 

111 IF (CG-6.0) 112*113*114 

112 PRINT 22-vl 

* 

22 OL FORMAT l 5^ *42H COVERGLASb I i-.COkRCCT-dUT WILL USE 6. 

113 MCG=1 

GO TO 2 - 0 2 

1 14 I F ( CG~9 • v ) 111,115*116 

115 NCG = 2 

GO TO ZOU2 

.116 IF(Cb“lS*C ) 111*117*119 * 

117 NCG=2 

GO TO 200- 


CHI LS) 


118. IF (CG-36.C) 111*119*120 



119 NCG = 4 

GO TO 2002 

120 IFCCG-54. 01111*121* 122 

121 NCG=5 

GO TO 2JCZ 

122 PRINT 2301 

2301 FORMAT (5X»40H COVFR-GLASS TOO THICK Vi ILL L'S" 54.CWILS) 
NCG = 5 

.GO TO 2002. 

2002 PH I £ 1 =0 • 0 
P H I E 2 = C • C 
P H IP 1 = 0 * 0 
PH I P2 = 0 *0 
PHIFP1=0*C 

. . .. PHirP2 = C.O 

PHIFA1=0*0 
P H I F A 2 = 0 • 0 
PHIE=0.0 
PH IP =0*0 
PHIFP= : *0 
P H I F A = . .C 
Pi -l I T = 0 • 0 
K-HCG 
J=NFS 
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HUM! I 


' )0 12 I = i *LEVF L£ 

I i ; = ( ; X - i ) L - V L ' L E + I 

Fill- } = { c ( I ) 0 A V S r ( I K ) ) + P H I E 1 

IJ={ J-l) 7 ' L F V F L : . + I 

c;UF? = ( 'r ( I )*LAY5*i;F ( I J ) >+PHl c :2 

123 continue 

equivalent electron flux- phie* 

Pi ; 1 1 " = PHIFI + r -i i F 2 

00 124 i = i»Lev;-LP 

1 k = ( , 2 - 1 ) *■ L E V t L P + 1 

r' i : 1 P 1 = ( r ( I ) 7 - u m Y ;*Ur ( I K j ) + P r i I P i 

I J = C J — 1 ) *■ L E V E L i -' + I 

Pi j I P2 = ( P ( i ) * 0 A Y A .0 P U J ) ) + Pitl P2 

124 CONTINUE 

1-r-EV £ v.iil VALENT PUvT^N FLUX -PHIP. 

PHIP= PHI P 1 + PHIP2 
;;0 125 I=l*LrVLFP 
IK = ( .<-! ) *LE VLFP+ I 

P H I F P 1 = ( FP( I )*LAY5*OP( IK) )+PHlFPl 
IJ=t J-1)*LEVLFP+I 
P H I F P 2 = t F P ( I )*0AYS*UP( I J) J+PHIFP2 

i 2 : 3 L0N1 I NUt 

1 - V ,E V lwoI VALENT 60 LAP FLAKE PROTON FLUX* 
PM I r P= PHIFP1 + P h I F P 2 
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DO 126 1 =1 * L E V L F A 

I ^ = ( - 1 ) «• L F V L r a + I 

Phlr Ai=(F/\( I ) '* L A Y 3 * ; > A l 1.0 )+PHlFAl 
I J= ( J-l ) *LEVLFa+I 

P H I F A 2 = ( f A t I )*DAY£*OA{ I J) > +P H I r A 2 
126 CONTINUE 

i-mev ' equivalent solas flare alpha particle flux. 

P H I F A = PH I FA 1 + PH I FA 2 

TOTAL i-MEV EQUIVALENT Jr ALL CHARGED PARTICLES. 

P n I T = PHIL + PH IP + PHIFP + PH I FA 
$ R I T E OUTPUT TAPE 3*20U4» PHIT* PHlt* PHIP* PHIFP, PH I FA 
RET URN 

20 J 4 FOR.; AT ( 26X *4lH TOTAL EOU I VALENT -1-iiEV FLUX=Ul0.3// 

1 2 X » 4 7 H FOU I VALENT 1 MEV ELECTRON FLUX=E 10 • 3 // 

222 X *47H EQUIVALENT 1 MlV PROTON FLUX = E 10 . 3 / / 

32*- X ? 47H EQUIVALENT 1 :*V SOLAR FLAPS PROTON FLUX = E10.3// 

4 2 0 X , 4 7 H 1 LEV SOLAR FLARE ’ALPHA PARTICLE FLUX*E10.3) 

l ; o r FOR:* AT ( 1HC , 441 1 GIVEN AS INPUT DATA 1 MEV EQUIVALENT FLUX = E 10 . 3 * / / 
1 /) 

2ND 

S' 'SPOUT ifi£ OGPADE ( PH I T t RHO * X I t V*MCELL *XPVtXP I ) 

0 1 HE MS IONVCC ( 101 ) tC I SC ( 101 ) ,VOC( 101) * A I SC ( 1 0 1 ) » X I (101) * V ( 1 0 1 ) 

DIMENSION IMAGE! 1717) » MSCALE(S) 

lF(\CcLL)2*2fl 
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1 CALL PlOTKNoCALEjIC * 1 C: * HI *10 ) 

NEAL INPUT TAP t 2 * 1 U C * NTuL* NPTS* LATE* TITLE 


PRIuT OUT OF 0 R I ‘ 

3 INAL LATA POINTS OF 

THE I — V Oj^VF. 

,-R I T F OUTPUT TAPE 

3 * 1 0 S * N T P L * N P T S * u A I : 

i T I TLE 

READ INPUT TAPE 2 

*101 * (Cl SC( I ) * V C C ( I ) 

, 1 = 1 »r<PT. c . ) 

Vi RITE OUTPUT TAPE 

3*1010 


WRITE OJTPJT TAPE 

3*102. (CISC ( I ) »VCC ( : 

I ) * I = 1 ...PTS ) 

THE ORIGINAL u;\TA 

P 0 I i T S A I 7 : E N 00 E X T K * 

: ■ P i j L A 1 r L • 

VOC ( 1 ) =C*U 



Aiscm=o.c 



DO 2 0; J = 2 * 1 a l 




A I 5 C ( J ) = 0 . 0 
VOCt J)=VCC<J-i ) +0.010 

* 

200 CONTINUE 
“1= 1 

DO 2 v 5 J = 1 *10 0 

IF (VCC ( I ) )2 02 9 *2031*2 C 2 9 
2C 29 IF ( VCC ( J ) -VCC ( I )) 2C2 * 203 *2030 
2F 3~ 1=1+1 

202 A I SC ( J ) = A I SC ( J- 1 ) + ( C I SC ( I ) -C I SC ( I - 1 ) )/<VCC< I)-VCC( 1-1} MM VOC ( J ) - 
1 VOC(J-l) ) 

•j? 0 TO 2 - 2 'j 
2031 A I SC ( 1 ) = C I SC ( 1 ) 

1=1 + 1 
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GO TO 2 0 2 0 

203 A I SC ( J)=CISC( I ) 

1 = 1 + 1 
2020 N* J 

I F ( C I SC ( I ) ) 204*2041 *205 

205 CONTINUE 

2 C 4 PZCtl = VOC ( N-l ) + A I SC ( N-l ) * ( VO C ( N ) -VOC ( M-l ) ) / C A I SC ( N-l ) -A I SC ( N ) ) 

GO TO 2: 42 

2041 VOC<N)=VCC< I ) 

Ai$C<N)=0.C 

PZCUsVOC (Ni 
N=N+1 

2042 OO 206 J=N*101 

A I SC ( J ) - A I SC (N-l )+ ( A 1 0 C ( N ) - A I S C ( N-l ) ) / ( VOC U ) -VCC ! N-l > ) * ( VLC (■-) 
l~VOC { N-I ) ) 

206 CONTINUE 

C THIS COMPLETES I-V CURVE EXTRAPOLATION. 

C CURVE IS NON DEGRADED. 10HM-CM :-.^T EE I AL + 10. CM' -CM MATERIAL. 

C 1 OHM— CM FIRST 

2 X M X = './ C C ( N P T S ) + • 0 5 
Y M X = C I 3 C ( 1 ) + • C 5 

CALL PLOT2 ( I.- AGE »X.iX * E.C * YMX * f;.0 ) 

CALL PLU T 3 ( lHu * VOC * A I SC > 1 0 1 ) 

X = 0 . J 
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y = o . ;• 

xi=; .. 

Y i = -■ • - 
X2 = 

Y2 = - 
X3=-.v 

Y 3 = v • - 
X 4 = _• • 4 

Y 4 = - • . 

M=~.w 

uR=2*: 

V = w • V 

VIiC = w • J 
VLOG = - • ^ 

Y L OG i - 0 • - 

vlog2=0.; 

I F{ RHO-1 0,0)207 *20 8 *225 
2,7 IF(PHIT - 1 • 0 F 1 8 ) 209 , 211 *210 
209 1 K ( P H 1 T - 4.CF12) 211*211*212 
2 1" R I - ( 2 • 3 - • 1 6 7 * LOG 1 C F ( PH I T ) ) 

CO TO 2217 
211 RI^l.O 

GO TO 2217 

2 12 1 =i • «-7 • 13/ ( ( i0*3 J « *3 • ji ) * t ( PHI T - 4. Oil 1 2 ) * *0 • 4 5 1 } 
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I 


I 


30 TO ^ 2 1 7 
io oh;;-c. 

208 IF ( PH IT- 1.CF15 >213*21^*214 

213 IF (PH IT-5* 2 FI2 )215*21?*216 

2 1/4 ? I - { 2 • 8 ' 6 -0*1323 * LOG 1 Or ( PH IT ) ) 

GO TO 220 1 
2 13 R I = 1 • 0 

GO TO 2221 

216 K 1 = 1* -2*78/( ( PH I T-5 • OE 12 ) **0* A 1 ) 

GO TO 2221 

VOLT j AT ISC IS r:0’,v DETER-UREu. 

1 OH-:-Ci : FIRST. 

221* iF(FnIT-3.0E14) 217*216*218 

2 17 I r ( P r! I T - 1*0 t 1 J ) 2 1 v * 2 1 v > 2 2 G 
v i p, \/ 1 c. r = ( .. , ; • 1 A- ; : 3 * L G C- 1 0 F ( P H I 7 ) - 0 • 1 ? 2 V ) 

GO To 223 
2 1 p V I ‘X = 0 * 0 
GO TO 223 

2 2 V U'lG - ( P H I T / 1 • 0 3 1 ^ J 

\/ 1 c = 0 8 * ( LOGO OF ( V LOG ) ) * * 1*8^) 

G G T v 2 2 8 

: ic riii-c:: 

2 2 21 I F ( 0 r i I T — .c • ( - ) aci *Z^2*2^^ 

2 2 1 I F { P i i 1 T -2 • 3 £ 1 2 ) 2 2 ^ * 2 2 . * 2 2 4 
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2 22 V IGC = ( •> * 2 1 2 9* LG51 JF ( PHI T )) -D. 1832 
GO TO 223 
223 V I SC = 0 • 0 
GO TO 223 

2 24 V LOG 1 = ( PHI T/2.5E13 ) 

v IoOl ;,.:G47# ( LOGIC F ( V’LCG 1 ) )**1 .h3) 

GO TO 223 

0 k f L A F 1 V L 0 K.E N- C I R Cu I T VOLT AGL o l GR A l/ / •• I i ON • 

t HE E 0 U A I I ON 0 U S E U F o R B u I h 1 CH< V . — C * . AND 10 0 Hn— C H RESISTIVITIES* 
223 IF (PH IT- i • 0 £ 1 4 ) 226*2 27*227 
226 IF ( °H I T -3.CE12) 228*226*229 
2 27 R V = t 1 • 7 7 9- 0 • C 3 8 8 * LOG 1 0 F ( PH I T i ) 

GO TO 2230 
228 R V = 1 • C 

GO TO 22^0 

2 29 VLOG2 = ( PHIT/3.0E12) 

RV= l 1.V-G.C222* « LOG10F < VL0G2 ) )**1 .67) 

GO TO 2230 

C I - v CURVES ARE NOW DEGRADED USING THE FOLLOWING COMPUTED FACTORS* 

C VALUE OF RADIATION RESISTANCE RR. 

2230 R R = V I SC/A I SC ( 1 ) 

00230 0=1*101 

V(J)=VGC(J)+ AISCtJ)*FR 
230 CONTINUE 
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C SHORT CIRCUIT CURRENT CORRECTION IS NOW MADE. 

DEL1AI «AISCU )#( 1.0-RI ) • 

IF (DELTA I -A I SC ( 1 ) 12301.2302.2 302 

2301 D0231 J=1.101 

XI ( J) = AISCU>- DEL T A I 
231 CONTINUE 

GO TO 23112 

2302 D023 1 11 J = 1 * 1 01 
XI ( J 1 =A I SC ( J 1 

23111 CONTINUE 

C DEGRADED SHORT CIRCUIT CURRENT. 

23112 DO 2311J-1.101 

IN — J 

IF { X I ( J ) ) 2332>2333»2311 
2311 CON T INUE 

C 

2? 3 2 DFL = ( X I U-l )*(VU)-V (J-l >"))/( XI ( J-l )-X II J) ) 

. If ?CD = V < J-1)+DL*L 
GO TO 2334 
2333 PZCD=V { J ) 

C 

2 3 34 SLOPE -(XI (N-D-XI (N-2) ) / ( V ( f> 1 ) -V C ; i-2 ) ) 
uO 2423 JaNtlOi 

xi c j ) = x i (f!-i) + {yu)-voi) ;*s Lppr 



i ’ »t*Nn in iHiiiiitirtti ii'i.imm i-Hi-Hic ihuoi iriiHH in. .< 


2423 CC':Tr;V . 

jLj 2 Pi;.. CIRCUIT VOLTAGE IS PZCU. 
iir.&X/.uf-'ti C !••>•?: Cl. ’CHIT VOLTAGE 15 PZCO*«V. 

v?c?- pzcw*hv 

2401 >FL T A V = PZ CD-VQCR 

rI>.L 1 ir 2A T f OF VOLTAGES IM THE ARRAY 
24.. 3 HO 2 3 0J=1 , 101 

V ( J ) -V ( J ) -DEL T AV 

25, v CO NT I is JE 

CALL PLOT 31 iHD,V*XI ,101) 

L033J=1 *101 
ir (V(J) ) J j » 32,^2 

32 X15C-XI ( J) 

GO ro 34 

33 COfiTliHlJF. 

34 .■ i * 1 

PwR= -1 ;0,o:).0 

2 3 30 PV,RT« V C f : i * X 1 1 .% » 

Ir (° "H-PWPI234.235.236 
236 N=N+1 

PwR= PwPT 
GO TO 2330 
23 5 VHP =V(N) 
x i mp =x I (/•;) 
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GO TO 45 

2 34 slope = ( x i (n-d-xi (nj )/<v(N)-v<n-d > 

237 VUP = V(N-1)+0*00015 

X 1 1 ■IP = XHN-l ) -SLOPE* (VMP -V(N-D) 

P'v.'R T = VMP *XIMP 

IF ( P ’ P T -PWR) 45*45*238 
2 ?« P‘\'R = PWRT- 

GO TO 237 

C FORMATS FOR PRINTING RFSUlTS. 

45 WRITE OUTPUT I APF 3*200" 

aRITE OUTPUT TAPE 3*2CC3* VOCR 
’.•/RITE OUTPUT T APF 3*2004* XISC 
WRITE OUTPUT TAPE 3*2031* XIMP 
WRITE OUTPUT T APF 3*2^02* V M P 
C MAX. POWER PF r R AT AT I ON • 

IF ( X P V ) 4001*4001*400" 

4'?o u; , p = x’' ■ * x p i 

DMP= t U ,4 P-V' T*xl f/ p ) / UO'P^IO^ ,0 
C Ol'ORT CIRCUIT CURRENT rjFG°.* OPTION 

1. IOC - ( CISC ( 1 )-XI r 0) /Cl SC( 1 } *ico,o 
"ITE OUTPUT TAPP 3 *2^-10* D”P *0150 

r i-oEouiru^i! percent cf maxim/*' power pfgrad/'-t icnf6. t/ 

1 i ' X * 4 5 h PERCENT ^F c -MORT CIRCUIT CURRENT OE 0 R A ;*•/• T I ON F 6 • 1 ) 
2. Cl FCR’-.AT ( 14H. ! *- f AX PUR F7.4) 
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2002 FORMAT { 1 4H V-MAX Pv/R F7.4) 

4001 DC 2312 J = 1 *10 3 * 2 

WRITE OUTPUT TAPF 3,2005* J*(XI(J)) 

2005 FORMAT (4H0 I(I3*?M) FT* 4) 

WRITE OUTPUT TAPE 3*2006* J*!V(JM 

2006 FORMAT <4H V(I3*8M) FT. 4) 

2312 CONTINUE 

1 2000 FORMAT (5 CMC IRRADIATION DEGRADATION OF SOLAR CELL I-V CURVE ) 

10 1C FORMAT (3ZHO ORIGINAL I -V CURVE DATA POINTS) 

2003 FORMAT (14H VOLTS OC F7.4) 

2004 FORMAT (14H AMPS SC F7 4) 

100 FORMAT ( I3*I2*5X*I6*4X*10A6) 

101 FORMAT! 4 (2F 10.0) ) 

102 FORMAT (lPHC CURRENT IN AHPS=F 10. 4/ 1 8K VOLTAGE IN VCLTS=FlQ.4/ ) 

. CALL FP LCT4 ( 23 *23H CURRENT) 

WRITE OUTPUT TAPF 3*104 
104 FORMAT ( 25X *?HV 6 L T S /// ) 

RETURN 

END 



DATA 


2.89F13 

3.1E2 

2.32E11 

1.91E6 

2.32E6 

0 • C E o 

1.91E6 

4 * 6 E 6 

X.09E5 

5 • 4E4 
. 1 • 23E3 
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2.35612 

7.78E11 

2.5E11 

1.8F2 

1 *0£ 1 

- 

3.56E8 

4.1E7 

2.I9E7 

1.23E6 

1.09E6 

1.09E6 

3.83E5 


• 

1.66E8 

4.1E8 

2.19E7 

1.23E6 

l."9E6 

1.Q9E6 

5.4E5 



8.2E4 

8.' :4 

8 • 2t 4 

2 • 7E^ 

. 3.2E4 

3 • 5E4 


- » < l3 


* 0 6 E 1 ^ 

4.09E8 

3.89F6 

1 • 2 3 E 7 

6.84E6 

8.47E6 

1.G9E6 

4.1E5 

3 • 12c5 

1 * 2 3 E . 

6.8E6 

5 • 4 E 5 

1.G9E6 

3.2E5 

6.3E5 

1.09F5 

5.4E4 

5*‘,t4 

1.36E4 

8. /oE5 

j * 4E4 


2.95E 4 

3.0E6 

2.46E6 

3.0E6 

4.9E6 

5 * 4 E 4 
6. 1E4 



26- . 

2 7 :• ( * 

2690 * 0 

2 50-0.0 

2500.0 

2500. 0 

25CO.0 

2500.0 




^ r 

A A 

A ^ Q 

A A 

2000.9 

23- . 

1 

• > 

■J • ' 





J 

3100.3 

2 6 0: . 0 

2700.0 

2600.0 

250-.' 

25CC .0 

2500.0 

2 5 w 

£ 0 J 9 ^ 

2300.0 

1 C u 0 * 0 

^ A 

0.0 

c.o 

0.0 


. *«* 

v- . 0 ! 

200.0 

1400.0 

2000.0 

2100. 0 

2100.0 

2 i - v • - 


2 0 Co *w 

20-0.0 

20C0.C 

1 C*oo .o 

C.O 

0.0 

1 • 

r 

^ w 

. -j 

.• . 0 

0.0 

0.0 

c.o 

0.0 


i Vrf . '-J 

i : . o 

14- 0 . 0 ■ 

1500.0 

1800. 0 

2 C 0 C • 0 

1000.0 



0 * 0 

^ 

O.o 

0.0 

C.O 

0.0 


O.o 

: *o 

O.o 

c.o 

"•C 

100*0 

1500.0 

<i ■ - - • • 

100 0 * c 







1 : - .. : . 

130.:.: 

14030.0 

137C0.0 

12500.0 

1150'-.. 0 

i 1000.0 

10400.0 

9 3 «. w • 

9500.0 

9 s O' o • 0 

9 5 0 r • " 

9 50° . 0 

c 4 0 - • 0 

9200. o 

8600.0 

7 \ ; . 

4 o : . . c 

* • c 

C. r - 

7000.0 

1100^.0 

12000.0 

11500.0 


I 0 4 0 0 • J 

9300*0 

95CG.0 

9500*0 

9500. 0 

9500.0 

9400.0 

9 2 0 2 • . • 

6 6 0 - : . -■ 

70 00*0 

4 0 C 0 . 0 

o.c 

n.o 

0.0 

o.o 


23.,.' 

o 2 ■ • • O' 

7200. 0 

79C0.C 

8000. 0 

7800.0 

7700.0 

76 uj.;j 

760 ... .u 

7 6 0 C *0 

7700.0 

f 3 0 0 • 0 

4C0w.C 

-.0 

o.c 

vJ • V 

^ 


— . 0 

0*0 

0.0 

C.O 

1700.0 

3 63--- .0 

Q : ^ ^ 

5700*0 

6300.0 

7000.0 

7400.0 

7000. 0 

4000.0 

; . • 

• ; 

C • 0 

n. r> 

c * 

0*0 

0.0 

0.0 

0.0 

■J « 

0.0 

C • 0 

'•J 9 0 

1600*0 

4000.0 

6400*0 

7100.0 

/ 0 -j ■j * c 

400 0*0 







1 43. 

0 6 . 

0 1C 

I-. 

. — ■ 

1365 

.0 0 

.4200 0. 

1300 
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,;0 113 Cl 1669 

HELICTFK 10 OH 

V CM CELL 

AT 30 DEG 

C 



'"■.1412 

0.0 • 1410 

o’.To 

0 • 140^ 

0.2000 

0.1398 

0.3000 

0*138 j 0 

.3pOo w • 1 3 & ~ 

0.3800 

0.1340 

0.4000 

0. ] 300 

0 . 42GC 

0.124C 0 

.4400 * j . 12 v 0 

0.450C 

0.1100 

0.4700 

0. 1000 

C.48C0 

0 • v 9 C C C 

. 4 90 0 C.CSuO 

0.5C00 

C.06CC 

0.5150 

0.C400 

0.5280 

0 • w 2 0 0 

C . 5 3 6 9. CO -• 0 

0.547 





U02 45.0 

18.0 1C.0 


365. 

q 

0.4200 

0.13C0 

003 15.0 

6.0 10.0 

2.0515 

1 


0.4200 

0.1300 


999 


EOF t SY 


1 CALL PL0TKNSCALE.10 »1C .10 .10 ) 

READ INPUT TAPE 2.1uC. NTBL > NPTS. DATE. TITLE 
C PRINT OUT OF ORIGINAL DATA POINTS OF THE I-V CURVE. 

WRITE OUTPUT TAPE 3 . ICO ,NTRL .NPTS .DATE .T I TLE 
READ INPUT TAPE 2.101 * (CISC(I) .VCC< 1 ) . l=l.NPT-U 
WRITE OUTPUT TAPE 3.1010 

WRITE OUTPUT TAPE 3.102* (CISC 1 1 ) *VCC( I ) *1 = 1 .NPTS) 

C THE ORIGINAL DATA POINTS ARE MOW EXTRAPOLATED. 

V0C(1)=C.C 
AISC(1)=0.0 
DO 200 J = 2 * iOl 
AISC(J)=G.G 
VOC( J)=VOC(J-l 1+0.010 

* 

200 CONTINUE . ; 

.‘•I- 1 

D0205 J=1 *100 

IF I VCC ( I > ) 2 C 2 9 *2031.2029 

2029 IF. (yO.Ct J)-VCC( I 1)202,203*2030 

2030 1=1+1 

202 AISCI JlfAISCf J-l) + tCISC( I l-CISCU-l 1 ) /{ vcc rn-vcc < 1-1 > )*<VGC( J) 
1 VQC<J-1 ) 1 
GO TO 2020 

2031 AI SC ( 1 1 =C 1SC< 1 1 

fir 


1=1+1 


00 TO 2020 


203 A. I SC < J ) *C I SC ( I ) 

1 * 1+1 
2020 N* J 

I F I C I S C I I ) ) 204.2041.205 

* 

.205 CONTINUE 

2C4 PZCU *VOC ( N-l ) +A I SC ( N-l )* ( VOC ( N ) -VOC IN-1 ) » / 1 AI SC IN-1 > -A ISC IN II 
GO TO 2-42 

2041 VOCIN ) *VCCI I ) 

AISC(N)=0.0 

PZCU=VOC IN ) 

N* N+ 1 

2042 DO 206 J=N»101 

AISC(J)=AISCIN-1)+ I A I SC I N I -A I SC I N- 1 ) ) / 1 VOC IN )-VOC (N-l ) »- 1 VGCIM 
1-VOC I N-l ) I 
206 CONTINUE 

C THIS COMPLETES I-V CURVE EXTRAPOLATION. 

.C CURVE IS NOW DEGRADED. 1GHM-CM MATERIAL + 1 0 . CHM-O: MATERIAL. 

C 1 OHM-CM FIRST 

2 XMX=VCC ( NPTS I+.05 
YMX=CISC(1)+.C3 

CALL PLOT 2 ( I, -AGE .XMX . 6*0 . YMX . 0.0 ) 

CALL PLOT 3 ( 1HU »VOC . AI SC . 101 ) 


X 2=^«0 


Y 2 * 0 • G 
X3«J.~ 

Y3*v- 
X4 = j • Q 
T4 = , 

RI=~.0 
kR = 0 • G 
RV=G.O 
V I SC = >- . J 
VLOG«~*£> 

VLOGl=0.w 

VLOG2=D.t 

IF( RHO-1C# 0)207 #208 #205 

207 IF(PHIT -1.CF15) 2C9#21G*210 

209 I F ( PH I T— 4.CE12) 211#211#212 

210 R I - ( 2 • 3 -G • 1 6 7* LOG 1 OF ( PH I T ) ) 

GO TO 2217 

211 \<l = 1.0 

GO TO 2217 

212 SI=l.~-7.13/( U0.CI**8«3l)*( (PHIT - 4.01112 ) **0.451 ) 


•30 TO 2217 
C 10 OHM - CM 

208 IF (PH IT- l.CE15>213.214,214 

213 IFIPHIT- 5. OE 12)215*215*216 

214 PI=(2.8C6 -0.1325 * LOGIOF ( PHI T ) ) 

GO TO 2221 

215 R I * 1 .0 

GO TO 2221 

216 R I * 1 • -2.76/J (10.T>**7.34)*< < PHI T-5 .0612 ) **0.41 > 

GO TO 2221 

C VOLTS AT ISC IS fiOW DETcRMINEG. 

C 1 CHM-Cr, FIRST. 

2217 1 F ( PHI T-3 .0E14 ) 217*216.213 

217 I F ( Prl I T- 1.0E13) 2iV.21V.220 

218 V I SC* ( » .01403*LQC10F ( PHlXI-G . 192y ) . . 

GO TO 225 

219 VUC=0.0 
GO TO 225 

22 VLOG= (PH IT/1.0E13) 

VI6C=(-.005*(LOG1 JF(VLOG) )**l.b5) 

GO TO 225 
C 1$ OHM-CK 

2221 I F CPHTT - 2 » ( >' 1 5 ) 221.2t2.222 
221 IF(Pl.ilT-2.5Eli >223.221.224 


If 


2 22 V !SC = ( - • 0 129# LOS i JF ( PH I T ) ) -0 • 1 3 32 
GO TU 223 

223 VIGC=C#C 
GO TO 223 

224 V LOGl = ( PH IT/2*5E1'3) 

V I5C=( - .0647* ( LOGIC F ( v'LGGl ) ) **1.43 ) 

GO TO 225 

0 k FLAT I V L 0 P.E N~C I RCU I T V 0 L T AG L 0 c. GH A i> / •. T i ON * 

L jAt'lE EOuAT IONS USED Fwk Lul h lUnt v > — C. , AND iu Chm-CM RESlSTIVlIIES* 

2 3 IF (PHIT- 1.CE14) 226 *227*227 
2 26 IF { °H I T - 3 • C L i 2 ) 222*223*229 
2 27 RV-{ 1. 779-C.C 588 *LCG10F( PH IT ) ) 

GO TO 2230 

228 RV=1.0 

GO TO 2230 

229 VLOG2= t PH I T/3 • 0£ 12 ) 

RV=( 1.«-G # C222*<LOG10F( VLOG2 ) )**1#67) 

GO TO 2230 

C i-v CURVES ARE MOW 0 FOR A DLL USING THE FOLLOWING COMPUTED FACTORS* 

C VALUE OF RADIATION RESISTANCE RR. 

2230 RR~ VISC/AISC ( I J 
uC2 30 J = 1 • 1 0 1 

V( J ) =VGC ( J ) + A I OC ( J) *PR 

230 CONTINUE 
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C SHORT CIRCUIT CURRENT CORRECTION IS NOW MADE. 

DELTAI »A I SC ( 1 ) * ( 1 •O-R I ) * 

I F (DELTAI -AISC(l) ) 2301 *2302 »2 3C2 

2301 D 0231 J=ltl01 

> 

XI<J)*AISC(J)- DELTAI 

231 CONTINUE 

GO TO 23112 

2302 D0231 11 J*1 t 101 
XI (J)*AISC( J) 

23111 CONTINUE 

C DEGRADED SHORT CIRCUIT CURRENT • 

23112 DO 2311J«1*101 
N* J 

IF ( XI ( J ) ) 2332 *2333 *2311 
2311 CONTINUE 

C 

23 32 D£L»(XZ( J-1)*(VCJ)-VCJ-1)) ) / C X I C J-l ) -X I ( J ) ) 
_P2CD»y( J-1)+D&L 
GO TO 2334 

2333 PZCD=V(J) 

C 

2334 SLOPE- (XI ( N-l ) -X I ( N- I ) / ( V (N-l I-VIN-2 ) ) 
y>0 2423 J S N 1 101 


XI ( J)»XI (N-1)«MV( J) -VI N-l ) i*SLOP£ 



2k2i CONTINUE 

UNDEGRADEO CPE.\ CIRCUIT VOLTAGE IS PECU. 
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GO TO 45 

234 SLOPE=(Xl (N-ll-XI (N) ) / ( V l N ) -V ( N-U ) 

237 VMP * VIN-ll+O. 00015 

X I HP « XI (N-l l-SLOPE*! VMP -VIN-ll) 

PWR T * VMP *XIMP 

IF(=WRT -PWR) 45.45*238 
23® P'.\'R = PWRT ■ 

GO TO 237 

C FORMATS FOR PRINTING RESULTS. 

45 WRITE OUTPUT TAPF 3.20C) 

WRITE OUTPUT TAPE 3.2CC3. VOCR 
WRITE OUTPUT TAPE 3.2004. XISC 
RITE OUTPUT TAPF 3. 2001. XIMP 
WRITE OUTPUT TAPE 3.20C2. VMP 
C MAX • POWER DFOF AT AT I CN» 

IF(XPV) 4001. 4001. 400" 

4 " 00 U:’P*XPV*XPI 

DMP* t U’lP-V? ‘P*X I HP ) /U V.o* 10" .0 
C SHORT CIRCUIT CURRENT DEGRADATION 

c i sc =(ciscd)-xircj /ci sci n*ioo.o 

"ITS OUTPUT TAPE 3. 20- lC.p'*P *01 SC 

23- 1” FORMAT (l"r,’7H PERCENT OF MAX I VU V POWER CEGRAD AT I ONF 6*1/ 
lli;X.45H PERCENT OF SHORT CIRCUIT CURRENT OEGRADATIONF6.il 
2 jC1 FCR 'AT(14h: I -MAX PWR F7.4) 


U 


' 


m 

*81 


2002 

4001 

2005 

2006 
2312 
20C0 
10 1C 

2003 

2004 
100 
101 
102 


104 


FORMAT < 14H V-MAX PWR F7.4) 

DO *312 J=1.1CI *2 

WRITE OUTPUT TAPE 3*2005» J.(XKJ)) 

F0RMAT(4H0 ) ( l 3 »PH ) F7.4) 

WRITE OUTPUT TAPE 3.20C6. J.tV(J)’ 

FORMAT (4H 7(1 3 *8H) F7.4) 

CONTINUE 

FORMAT ( 5CH0 IRRADIA T ION DEGRADATION OF SOLAR CELL I-V CURVE ) 
FORMAT ( 3ZHO ORIGINAL I-V CURVE DATA POINTS) 

FORMAT (14H VOLTS GC F7.4) 

FORMAT (14H AMPS SC F7.4) 

FORMAT! 1 3.12 »5X. 16. 4X.10A6) 

FORMAT (4 (2F10.0I) 

FORMAT ( 18HC CURRENT IN AMPS=F1C.4/18H VOLTAGE IF V0LTS=F1 0.4// ) 
CALL FPLC T4 ( 23 .23H CURRENT ) 

WRITE OUTPUT TAPF 3.104 
FORMAT ( 25X.9HV 6 L T S /// ) 

RETURN 

END 
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APPENDIX B 


DATA 

on 

2.89F13 2.35E12 7.78E11 2.5E11 1.06E10 4.09E8 3.69E6 2.95E4 

3.1E2 1.8F2 l.OEI 

oi8 ; 

2.32E11 3.56E8 4.1E7 2.19E7 1.23E7 6.34E6 5.47E6 3.0E6 



1.91E6 

2.32E6 

O.OEC 
1.91E6 
4*6E - 

1.09E5 

5.4E4 


a • 2 3 E 3 


1.23E6 

3.83E5 

1 .66E8 
1 .23E6 
5.4E5 

8.2E4 

2.7E4 

2.7E3 


1.09E6 

4 • 1E8 
1.09E6 

8.2E4 

3.2F4 


1.09E6 

2.wl/ 

1.09E6 

V 

8 .2E4 
3.5E4 


1.09E6 

1.23E7 
1 • 09F6 

1.09E5 

1.36E4 


4. c 5 


6.8E6 

8.2E3 

5 .4E4 
8.76E5 


3.12l3 

5 • 4E6 
6.8E5 

5 • 4E*» 
5 • 4E4 


2 • 46t6 

3.0E6 

4.9E6 

5.4E4 
6 • 1 E 4 


O.G 0.03 0.13 0.30 1 • ’ 5 2. 70 4.13 5.30 

6.15 7.30 7.8C 0.0 0.02 0.08 0.20 1.02 

2.50 3.92 5.15 6.10 7.30 7.80 0.0 C.O 

-.03 0.10 C. 75 2.05 3.36 4.60 5.70' 6.80 

7.70 C.C O.' 0.0 C.02 0.47 1.55 2.83 

4.13 5.3C 6 ”0 7.5C 0.0 0.0 0.0 0.0 

-#•25 1.10 2rl5 3.30 4.60 5.85 7.0 

0.0 0.0 0.0 C.O 3000.0 370^.0 3-30.0 3100.0 
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2500.0 

25C0 .0 

2500. 0 

25CO.O 

2500.0 

28 — . 7 

2 70 C . ■> 

260C .0 







r** t i 

0.0 

O.n 

o.O 

0.0 

2000.0 

23 -* r • 

loo . j 
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